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The best-known theoretical works on boundary-layer problems, especially in the case 
of compressible flow without or with heat transfer, are related to the laminar boundary 
layer, although the turbulent boundary layer is, in practice, often more interesting. The 
laminar boundary layer is more easily accessible to theoretical treatment because clearly 
defined relations exist between the viscosity m and the shear stresses r. In the turbulent 
case, empirical relations must be introduced. Therefore, attempts to get exact solutions 
are not worthwhile, while efforts to obtain approximate solutions, based for instance on 
the von Karman-Pohlhausen principle (Z. Angew. Math. u. Mech. 1, 233 and 252, 1921) 
of utilizing integral conditions, appear to be appropriate to this problem. In the last few 
years the accuracy of such approximate solutions for the incompressible case was notice- 
ably improved by the application of a new energy integral condition in connection with 
a new empirical law for the dissipation in turbulent boundarv lavers. stated by J. Ilotta 
(Ingr-Arch. 20, 195, 1952) and E. Truckenbrodt (Ingr-Arch. 20, 212, 1952). The empirical 
laws for dissipation and for turbulent wall friction, which are needed in the present approxi- 
mate theory, are formulated on the basis of available measurements for incompressible flow. 
Generalization to the compressible flow with heat transfer is made from physical considera- 
tions. Calculated results agree satisfactorily with available experimental data. Some pos- 
sibilities for improving as well as for simplifying the approximation theory are outlined. 

1. Introduction 

In principle, problems of compressible laminar boundarv layers may be solved by using 
the partial differential equations for the velocity and temperature 1 fields, and the thermodynamic 
relations between viscosity /x, density p, pressure p, and temperature T. Exact solutions for 
special cases have been worked out by several authors [5, 6, 7, 8]. 2 Nevertheless, approximate 
theories of the von Karman-Pohlhausen type have been successful because of their easy adapt- 
ability to the most general case of flows with arbitrary pressure gradient in flow direction 
[9 through 16]. 

An example is the case of a eompressible turbulent boundary layer, which is more important 
in practice because the controlling viscosity is no longer a material property, such as the molec- 
ular viscosity, \x\ rather, it is like so-called "apparent viscosity/' fx a , which is essentially 
determined by the empirical laws of turbulent mixing motions. In view of the empirical 
character of these laws, attempts of exact solutions of Prandtl's boundary layer equations 
are indeed inappropriate. Hence, approximation theories of the von Karman-Pohlhausen 
type are indeed the only ones which are well adapted to the problem of calculating turbulent 
boundary layers. 

The first solution was given by E. Gruschwitz [17] in 1931. He introduced a one-parameter 
family of empirical turbulent velocity profiles with parameter H, and, as the second param- 
eter of the problem, the boundary-layer m.omentum-loss thickness <5 2 . Gruschwitz then 
derived an empirical energy law. This was used together with the known von Karma n integral 
condition for the momentum within the boundary layer, to determine the two parameters 
H and d 2 as functions of x for a given pressure gradient. 

The so-called Gruschwitz energy law could not be verified in a general manner and appears 
therefore somewhat doubtful. Moreover, it is restricted to the incompressible case. A definite 
improvement of the approximate theory for turbulent boundary layers was accomplished 

i This work was supported by the U.S. Air Force, through the Office of Scientific Research of the Air Research and Development Command. 
2 Figures in brackets indicate the literature references at the end of this paper. 
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when several authors [9 through 15] independently (between 1935 and 1946) showed that 
an integral condition of the energy within the boundary layer may be derived from the Prandtl 
boundary-la^^er equation in the similar way as for the momentum. Some years later (1950) 
J. Rotta [3] derived from measurements the empirical law for the dissipation within the turbulent 
boundary layer, which was needed for the applications of the new integral energy condition. 
At about the same time the empirical law for the turbulent wall friction occurring in the mo- 
mentum condition was improved by new accurate measurements of Ludwieg and Tillmann 
[18], such that it may be generally valid for positive and negative pressure gradients in flow 
direction. 

Now, on the basis of the momentum and energy integral conditions on the one hand 
and the two empirical laws for the turbulent wall friction and dissipation on the other hand, 
a two-variable method approximate theory of boundary layer could be developed, similar 
to that formerly worked out by several authors for laminar boundary layers [9 through 16]. 

For the incompressible case, the integration method of E. Truckenbrodt [4] (1952) (similar 
to that developed by A. Walz [19] in 1941 for the laminar case) appears to be the simplest. 
However, applications to compressible turbulent boundary layer failed, because the semi- 
empirical laws for wall friction and dissipation could not be set up, since before that time 
measurements of turbulent boundary layers at high Mach numbers were not available in 
sufficient generality. 

Basing his analysis upon simple physical facts about the turbulence mechanism the 
author was able to derive a generalization of the empirical laws for wall friction and dissipation 
to include compressibility, [26, 27]. 

The inclusion of the case of heat transfer, however, implies a new difficulty insofar as the 
integral energy condition derived from Prandtl's equation considers only an energy flux along 
the body surface but does not relate to an energy flux through the wall. In the latter general 
case we must start from the general differential equation for the temperature field. In addition 
it must be checked whether the above compressible considerations would include heat transfer 
at the wall. 

As suggested by L. E. Kalikhman [13] for the case of laminar boundary layers the partial 
differential equation for the temperature field can be transformed by partial integrations into 
an infinite system of ordinary differential equations similar to that obtained by partial integra- 
tion of PrandtPs boundary layer equation. 

Hence, two infinite systems of ordinary differential equations with a twice infinite number of 
unknowns are available which appear to be completely equivalent to the two partial differential 
equations for the velocity and temperature field. 

Some general properties of these systems, which had not yet been outlined in the literature, 
will be discussed. 

As a first step of evaluation of the possibilities involved in these systems of equations, a 
"two-equation, two-unknown" method is investigated. This is about as simple as the known 
methods for the incompressible case [4, 19]. Nevertheless results obtained are of sufficient 
accuracy for engineering purposes. 

The application of the new method is illustrated in some examples. Directions for improve- 
ment as well as for further simplifications are indicated. 

2. Partial Integration of Prandtl's Boundary-Layer Equation for a Compressible 
Two-Dimensional Steady State Flow Along an Insulated Wall 

2.1. Derivation of an Infinite System of Ordinary Differential Equations 

As proposed by Leibenson [9], Golubew [10], Sutton [11], K. Weighardt [12], and E. Truck- 
enbrodt [4], PrandtPs system of partial differential equations for a compressible two-dimen- 
sional flow along' an insulated wall may be transformed into an infinite system of ordinary 
differentia] equations by partial integrations. We will give here a new derivation of this 
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system which is appropriate for our purposes. 

Let x 9 y be the coordinates, the aj-direction being taken everywhere tangent to the wall. 
Further u and r are the components of the velocity vector, respecting in the ./; and y directions, 
m the viscosity coefficient, 6 the density coefficient. Then PrandtPs boundary-layer equations 
may be written as: 

d( —} 

du . du du& . \ dy/ /1N 

pu vx +pv *r pm te + — *f-' (l) 

dfai) d^) =0 (confcinuit equation) . (2 ) 

Of oy 

Here the subscript 5 denotes that the quantity is referred to the edge of the boundary layer. 
We note that this system, is valid even for turbulent boundary layers, if the velocities u and v 
as well as p are considered to be the mean values over time and if p is the effective value p = p e , 
which consists of the molecular viscosity /^, the apparent viscosity id,,- Thus, 

/x e = M + M« • (3) 

For simplicity, however, we prefer to use in hot h cases, laminar and turbulent, the same symbols 
u, V, p, and ji, knowing that in the turbulent case we have p = p € - 

First we multiply eq (1) by u" and add to it eq (2) multiplied by u v+1 /(v-{-l). 

Partial integration from y=0 to y= oo then yields after elementary transformations, 

' (pM +2 f,)+g,'?^+e,=0 (F=0,l,2,... f »).» (4) 



v + 2 1 \HQ^5 JV) I \v 

PbU\ ' 2 ax 
Herein the following notations are used 



J \u s / du\p 6 ui/ 

*-/."D£-s(s)1* (6 > 






(8) 



From the structure of the integral terms (5) through (7) we see that for practical purposes it 
will not be necessary to extend the integration process beyond a certain distance y = 8, for 
which we have u=u 5l p = ps, r = 0, because for y^>8 all integrands became practically zero. 
The layer thickness 8 is somewhat vague in view of the asymptotic behavior of u(y) at the 
outer edge. For approximate theories, however, it is very useful and may be defined as wall 
distance y at which u reaches the value % m with an accuracy better than 99 percent. Hence 
we put 

0.99tt„<^<>co. ( 9 ) 

2.2. Discussion of the Infinite System of Equations 

Putting v=0 we obtain the usual momentum integral condition of von Karman-Pohlhausen, 
for v=\ we have Weighardt's energy integral condition. Only these two cases have been 
considered up to the present, but it is clear that all other remaining integral conditions of the 



3 v may have also values between — oo and + 00 , even nonintegers as 0.1, 0.2, . . . and so on. 
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infinite system (4), the so-called higher momentum conditions (v=2, 3, 4 7 . . . ,00), would be 
important too, since only the complete system is equivalent to PrandtFs equations. 

Let us first investigate the meaning of the equations of order v, higher than unity. 

First eq (4) is rewritten in a more convenient form. At the outer edge y = 8 of the boundary 
layer, we may assume adiabatic conditions of state. According to the basic laws of gas dy- 
namics we have between the density of p&, the velocity Us, and the Macli number Ms=us/as 
(a 5 =local sound velocity) the following relations (see, for instance, It. Sauer [20]): 

lf = ip m . (10) 

p8 dx ud dx 6 
Substitution of (10) into (4) gives 

r,+^/„[2+F-2KJ4^]+*=0 (11) 

where the prime denotes differentiation with respect to x. 

The integral terms e vy f P , g v represent the unknowns of the system. They contain the 
velocity and density distributions u(x,y)/us(x) and p(x,y)/p 8 (x), respectively, which are indeed 
the real unknowns of the problem. 

The only way to obtain a solution of the problem was to introduce approximations for 
u/us and p/p 5 according to von Karman-Pohlhausen's idea, for instance by means of a one- 
parameter family of experimental curves. In the incompressible case (p/pa=l) the problem 
was indeed solved in this way by taking two equations from system (11) (for v=0 and v=l)f 
using as unknowns the boundary layer thickness 8 or one of the functions e v ,J v , g„, and, as the 
second unknown, a shape parameter H of the velocity profile. In this case the ratios e v /5, 
fv/5, gv/5, as well as e v /f v and g v [f v , are only functions of the shape parameter H. 

The compressible case involves new variables p and T which are related by 

-*-4 (13) 

Pb 1 

since dpldy = 0. Additional equations will be needed to make the problem determinate. 
This is done for the thermal boundary layer. While we will return to this problem in section 
3 of this report, let us first discuss the principal properties of the infinite system (11). 

The order number v may for instance be increased toward °o . In this case we have from 
equations (6) and (7), 

[f v (v+2-M$)+g v U m =fM-M$), (14) 

and 



f 

J oc 



because 



-T^dy, (15) 

Jo Psu s y ' 



(*)"*-(»)' =0, for0<^<l, (16) 



< The original Pohlhausen method only used the equation ^=0 (momentum equation) from this system. As the second equation, Pohlhausen 
introduced the so-called wall condition, 



(*w),.r£- (12) 



which follows directly from Prandtl's eq (1) for the point y=0. It is obvious that this condition, valid for the wall point, would be doubtful if used 
for fixing the shape parameter. Indeed essential improvements of results were observed when this wall condition was neglected and replaced by 
the equation for v=l (see, for instance, [16]). 
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whatever the functions u/us and p/p 8 may be in dot ail. 

The integral expression (e v ) v ^ 00 is found to be zero as follows: 



it 
For 0<— <1 it yields 



[ (i+ ">©1-,.-°- 



Since furthermore the shear stress r and its derivative dr/d?/ are vanishing at y = 8 (n=u 8 ) for 
arbitrary velocity profile u/u 8 , we have 

(e„),->«,-»0. (17) 

Hence in the limiting case of y = oo, system (11) reduces to 

/:+/„(l-M?)5=0. (18) 

This equation may be called the mass-flow equation and is nothing less than an averaged form 
of the continuity equation. With the aid of (10) it can immediately be integrated to 

/coP«w«= (8 — 5 1 )p 8 u 8 =const, (19) 

where 



Jo \ piUj '' 



(20) 



is known as the displacement thickness. This relation has not previously been used in approxi- 
mate boundary layer theories. It expresses that the mass flow between the stream lines y = 8 
and y=^di must be constant. We have not ascertained whether relation (19) will be very 
useful in practice. But it is clear that if the velocity and density profiles are once determined 
in the inner part of the boundary layer (thereby determining the value 5^, the total boundary 
layer thickness 8 is determined also by relation (19). 

We have found that if the two-equation system v = 0, v — l with two unknowns is used, this 
relation is mostly violated. Nevertheless the essential properties of the boundary layer, for 
example, wall friction and separation behavior, appear to be correctly described. 

From this consideration about the mass-flow equation as limiting case for v = co, we may 
learn that the higher the order v, the less details of the velocity profile will be included in the 
corresponding equation, because the term e v , which contains the shearing stress (and with 
this the entire boundary layer effect), approaches zero, with v— >oo . 

It is therefore not surprising that equations corresponding to low order, for example the 
equation for *> = and v=l f still give a very good description of the properties of the boundary 
layer if a rather physically reasonable family of velocity profiles is introduced into the equations. 
A good way in principle as suggested for the incompressible laminar boundary layers a long time 
ago by several authors [9,10,11,12], would be to restrict to a finite number of equations of the 
system, preferably to those of low order, say up to v — 5 or 10, and to introduce a polynomial 
expression for the velocity profile 



u 8 



^[^ly^W^Cx),^^)...., ff,(*)] (21) 



with v coefficients H. Then v-\-\ ordinary differential equations will allow the determination 
of the v coefficients H v and 8 2 as functions of x. In this case no detailed assumption about the 
velocity profile, except some trivial boundary conditions, must be made. 

It appears worthwhile to study this problem with new efforts, because its solution will 
be to a large extent equivalent to the solution of PrardtFs boundary layer equation, but needs 
no restrictions as to the prescribed pressure distribution p(x), and will in principle be applicable 
to the more important case of turbulent boundary layers too. Since we want to include the 
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compressible case with heat transfer at the wall, we must simultaneously consider the temper- 
ature boundary layer. 

The infinite system of equations derived from PrandtPs equation contains no terms which 
take into account an energy flux through the wall. Therefore, only the momentum equation 
p=:0 and the mass-flow equation v=°° of this system can be used in connection with heat 
transfer problems. All other equations are worthless in this respect. The essential step to 
get out of this difficulty is to include the partial differential equation for the temperature field 
in this problem by transforming it into an infinite system of ordinary differential equations 
similar to that obtained from PrandtPs equation. 

3. Derivation of Integral Conditions From the Partial Differential Equation for 

the Temperature Field 

The general partial differential equation for the temperature field in a two-dimensional 
steady state flow may be written as follows, if the enthalpy i=c p T is introduced: 

di bi dp . /du\ 2 d /A di\ , OON 

Here X denotes the heat conductivity (kcal/m sec °K),c p the specific heat at constant pressure, 
u, v must satisfy the continuity equation 

^) + ^) ==0 . (23) 

dx dy 

Now, if we multiply eq (22) by i, and eq (23) by z' v+1 /V + l and add the multiplied equations, 
we obtain 

5( P m" +1 ) . d(pvv +1 ) , . ^ f du 5 . . fdu\ 2 . v d /X bi\. \ , OA . 

-^T-+-^~= (v+1) { - Up ^ lx~ lV +» h*y) l+ Zy fe Zy} 7 ' (24) 

By partial integration of eq (24) with respect to y from y=0 to y= CG and by substitution for 
v from eq (23) we have 



= (jr+l) I -p 8 us -^J o (uv-u s %)dy-p s us j£ udl*+j o i v d [- X - ^]+ J Q M* (If) <fy J- ■ (25) 

The terms with the factor 5 vanish because of the following known relation in gas dynamics • 

1 di s v 2 8 1 du& , , 

is dx i 8 u s dx ' 

which follows after a short intermediate calculation from 



^I=(jt-1)M|, k=- p -> (27) 

is 



U _A— (h 1 \ 71 /f 2 h — °P 



and 



9 Vf* 2 ?/, 

M$=i±7 ~ -t^T ; M%=-*) a*=critical sound velocity (28) 
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d\I b 1 <hu ( t-1 . 



■sSO+t 1 *) 



M 5 dz 
Using notations for the integral terms: 

J o P5^5 L\ P / J ' y ~ J (» Pa'W« L \ P / J V 
(oo can be replaced by 5, because the integrand is almost vanishing for 2/^5) 

Jo\P/ rw \uj 

J Vp / b^L^ 5(?//^) J 

dw r dul fj,c p D -d ul . 

w = yu 3— j ~—^=t J r =z constant rrandtl number 

77^=-=— =( because ^-=0 h 
equation (25) may be rewritten as follows: 



with 



and 



T=MT~' Tu 



d<t> v . , du s L-M% 

-1 — r<Pv -1— 

ax ax u 8 



.(„+!) *-lMf : 



pmV 



p*-i 



(29) 

30 

(31) 
(32) 

(33) 
(34) 

(35) 



</>,., €„ iiikI y„ (with v as an arbitrary Dumber) arc an infinite number of unknowns in I his infinite 
system of ordinary differential equations. The real unknowns, however, arc the distributions 
u(x,y)/us(x) and T(x,y)/Ts(x) which we are trying to determine by moans of the infinite number 
of integral conditions as expressed in the system (35). Now, before going to the solution, let 
us transform it into a somewhat more convenient form. First a new variable is introduced: 



<k 



where 5 2 is the displacement thickness, defined by 



JoP&u>s\ u s / 



(36) 



(37) 



It is to be remembered that t\ v is finite when M 5 -^0, while 4> v has the disadvantage of becoming 
zero when M 5 -^0. Furthermore, 



With those notations (35) can be rewritten as follows: 






(,=0,1,2,. . .,»). 



(38) 
(39) 

(40) 
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In principle a possible way to solve the problem of determining the distribution u(x,y)/u?,(x) 
and p(x,y)/p8(x) is to use polynomial expressions which fulfill the boundary conditions of the 
velocity and thermal boundary layer, by choosing j coefficients which satisfy j equations of 
the system (40). It is obvious that in the present case the procedure is more complicated 
than in the case of system (11). However, it may be possible to overcome by taking the solu- 
tion from the exact theory for the temperature profile T/Ts or the density profile p/pb=T§/T. 
Although the solution is exactly valid only for the special case of the flat plate (dp/dx = 0), it 
may be considered as a good approximation even for the more general case with dp/dx^Q 
(see for instance H. Schlichting [21, p. 282]). An essential property of this solution for T/T d is 
that it depends not upon x but only upon u/u§ with the Mach numbers M§ and the heat transfer 
coefficient c T as parameters. 

With reference to a paper of von Driest [23] this solution for T/T 8 for the case dp/dx=0 
can be written as 

Y^+cAfi-^+r — MUl-h) (41) 



irith 



Kl' P ')~l' *($ P >)~($ *»0*<P,<1, <«> 

J 7 rp 

--— \ — ^=recovery factor, (43) 

J- o — J 8 



r 

'-0—J-8 



T r —T w 

Cr—^-Tp — -=heat transfer parameter. (44) 



Thus the temperature and density profiles occurring in the integral conditions (36) to (39) 
are uniquely correlated to the velocity profiles u/u 8 , which now are the only unknown functions 
in the problem. So far all considerations and derivations are valid both for the laminar and 
turbulent case. For the solution of the laminar case, we have only to introduce an expression 
for the velocity profiles. The viscosity, as principal painmeter of all boundary layer theories is 
known to be dependent only upon the temperature T, for instance according to Sutherland's 
law. 

mi T+C\Tj {45) 

(with C'=113°K for air, T x being a reference temperature) or to its more convenient approxi- 
mation 



HtJ 



mi 
with 

Hhi+eyF, (47) 

(co~0.80 for temperatures T between 233° and 373°K, co = 0.5 for very high temperatures, 
a)=1.5 for very small temperatures near the zero point). 

In the case of turbulent boundary layers which are mainly considered in this report, an 
empirical relation for the effective viscosit}^ (as the sum of the apparent and molecular viscosity) 
must be introduced besides the assumption for the velocity profiles. Before this is done, a 
simplification of the systems (11) and (40) will now be considered. 
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4. Special Case of a Two-Equation System 

4.1 General Considerations 

As a first stop of evaluation of the possibilities involved in systems (11) and (40) we will 
make a so-called one-parametric representation for u/us, viz: 



u u ( y tts \\ 

u 8 u 8 \8(x) V 



(48) 



where H(z) is the form parameter of the velocity profile (as one unknown of (he problem) and 
d(x) is the boundary layer thickness (as second unkown). Therefore, to he correct, t lie relations 
(48) should be denoted as a two-parametric one. 

It must be pointed out that 8(x) may be replaced by another characteristic length of (he 
boundary layer, preferably by the momentum-loss thickness 8 2 (;x), thus avoiding any assumption 
as to the point where u/u 8 reaches the value 1. 

Since we have now decided to use only two unknowns, II (x) and 8 2 (x), we need only two 
equations besides the empirical law for effective viscosity. 

We choose the simplest equations (namely those for v = 0) from both the system (11) based 
upon Prandtl's equation (1) and the 4 system (40) based upon the general temperature field 
equation (2). 5 

These equations are: 

a. The momentum integral condition from system (11). This is the only equation of the 
system (11) which, besides the mass-flow equation, is valid even in the case with heat transfer at 
the wall. 



/;+/„ g+2-A/f]5 + , =0 



with 



pM 



JoPsUiL 1l S J J 2 ' 

Jo\ PdUs/ 



(49) 

(50) 
(51) 

(52) 



b. The energy integral condition from system (40) 



77 +77o$ — +y-=0 

U 8 2 



with 



T?0 = 



1 



Jo Us\ pj 



(k-l)M 2 8 



$=l+(fc+l)Mf -^ (independent of v), 

0-2 

PiU's \J« t, c \u s / 2s J 



(53) 

(54) 
(55) 
(56) 



5 The value of the simple mass-flow equation (19), which may also be obtained from system (40) by putting *>= - «> , appears to be limited, but 
has not yet been studied in detail. 
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CT ±T J- W 



(57) 



For the practical application of the two systems (49) and (53) it was found to be very 
helpful to evaluate once for all the universal functions 8 U 8 2) 770, $, and ^ for the interesting- 
range of the parameters H, Ms, and 6 (or c T ). Such an evaluation will be based upon a certain 
physically reasonable assumption about the velocity profile, as represented by (48). These 
assumptions for u/u 8 must of course be different in the laminar and turbulent cases. In the 

turbulent case, the empirical laws for the wall friction r w and the dissipation integral rda 

must also be introduced. The evaluation for the laminar boundary layer contains no diffi- 
culties, and will be dealt with in a separate report. Only the more complicated case of the 
turbulent boundary layer, which needs some special considerations will be treated in this report. 

4.2. Two-Equation Method for Turbulent Boundary Layers 

We now consider: a) A suitable expression for empirical turbulent velocity profiles, 
b) an empirical law for the wall friction in a compressible turbulent boundary layer with heat 
transfer, c) an empirical law for the dissipation in a compressible turbulent boundary layer 
with heat transfer, and d) a method of solution of the system of equations (49) and (53). 

a. The Turbulent Velocity Profile 

Detailed investigations of many authors [23, 24, 31 have shown that the turbulent boundary 
layer can be divided into several (usually into two) regions with respect to the direction normal 
to wall, in each one of which different laws for the influence of the outer pressure distribution 
and the Reynolds number are valid. Thus a general representation for the velocity profile 
requires a somewhat complicated formula. If onl3 r mean properties (integral quantities) of 
the boundary layer are considered and put into correlation as in our case, it will be sufficient to 
have an expression which contains less details but includes the whole range of the boundary. 

Such a formula, which is proved valid even at high Mach numbers, is 



Ui 5 

where k is the form parameter with values approximately 0.1 and 0.7, and is a unique function 
of any H to be chosen. 

A disadvantage of (59) is of course that at the wall y = the value (du/dy) y= Q is infinite for 
all 0<&<O and that furthermore, at y/8 = l, dujdy does not vanish as is physically required. 
But we have shown in previous publication [26, 27] that these facts do not essentially influence 
the integral terms of the approximate theory, since the wall-friction law and the relation between 
the velocity profile u/u& and the shear stress profile t/t w are separately introduced as semi- 
empirical laws. 

As an example of the suitability of this power law even for Mach numbers up to about 7 
and with heat transfer at the wall, we refer to figure 1 which is drawn from experiments of E. E. 
Wilson [25]. This simple power law is very adequate to the purpose of this approximate theory 
with regard to the possibility of closed evaluation of the integral terms. Figure 2 shows the 
order of magnitude of the errors in the ratio (d 3 /d 2 )i plotted against (di/d 2 )*, 



(<«-r5[-©» 
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Figure 1. Measured turbulent velocity profiles at a Mack num- 
ber M d = 6.8 and cooled wall (- = 0.840 (X) and 0.525 (#) X 

from R. E. Wilson [25]. Figure 2. Ratio (8 3 l8 2 ) = H plotted against 

the ratio (5i/5 2 ) t for turbulent velocity 
profiles. 

Comparison of experiments (according to J. Rotta and 
K. Weighardt) with approximation power law — = 

(f)«H). 

caused by the approximation of the experimental velocity profiles with the power law calculated 
for the incompressible case. We see thai the maximum error is about 1 percent. The neglect 
of a certain influence of the Reynolds number on the velocity profile causes errors less than 1 
percent. Therefore, the use of the power law (59) for our approximate theory appears to be 
justified. 

Some improvement at the outer edge region of the boundary layer can be reached without 
much complication by taking the following form: 



it 
u 8 



=tanh [b log (1+ay)], 



(GO) 



where a and b are two form-parameters, comparable to II X and IL. a may be used for the pur- 
poses of the adjustment of scales in the one-parametric theory, and the adjustment of wall 
friction in the two-parametric theory. The relation (60) avoids the use of 3 which is physically 
vague, and yields automatically an asymptotic behavior near the edge of the boundary layer. 
At the same time it prevents an infinite slope du/dy at y = 0. 

Investigations of the practical importance of expressions (60) (which of course, excludes 
the possibility of closed integrations) are in progress. 

b. The Wall Friction Law 

For this we may refer to the author's publications dealing with the case of compressible 
turbulent boundary layers without heat transfer [26, 27]. We now have only to investigate to 
what extent the relations found there will be influenced by heal transfer at the wall. 

As a summary of these earlier investigations we can state that the laws for turbulent mixing 
motion found for the incompressible case, are valid in the compressible case, provided that the 
fluctuation velocity is less than the sound velocity, in other words, if the velocity fluctuations 
are not accompanied by density fluctuations. Now, since the turbulent fluctuation velocity 
in boundary layers is found to be small, say 10 percent or 20 percent of the main velocity u& 



508986—50 5 



63 



0/. 

at the outer edge of the boundary layer, M 8 =— may be increased to about 10 or 5, respectively 

a 8 i 

until the fluctation motion becomes compressible in character. Therefore we conclude that up 
to about Mach number M=10 or 5, respectively, the laws for incompressible fluctuation may 
be considered as valid to a first approximation. Only the flow medium properties /x and p are 
changed. These, however, are taken into account by a suitably chosen Reynolds number. 

As to the effect of heat transfer at the wall, this will change only the flow medium properties 
and therefore can be included in the compressibility effect and in the Reynolds number, which 
itself plays a decisive role in the wall friction and dissipation law. 6 

The momentum law (V=0 from system (11)) contains the term r w \p 8 u\ which is the 
ratio of the wall friction t w to twice the stagnation pressure. This ratio can be rewritten 
as follows: 

/du\ r d(u/u 8 ) l r b(u/Ui) 1 

_^ J\*y) w ± dyft, lj^W^)} w ^ {HM (^f^l-^V (61) 

p 8 u\ p b u\ p 8 v 8 8 2 R 8 (5 2 )i JoU 8 \ uj 



A definition of a Reynolds number 



R PJ^l (62) 

Mm 



is suggested which contains the viscosity p. w at the wall and, by means of the integral term 
6 2 , an averaged inertia term p 8 u 8 8 2 . Indeed the momentum integral condition (49) brings into 
relationship the wall friction and the averaged inertia forces. The pressure forces are not 
of interest in this connection. Any other definition of the Reynolds number appears to be not 
suitable in the present problem. In the incompressible case these detailed considerations are 
of course not necessary. 

In the numerator of eq (61) we have the dimensionless slope of u(y) at y=0. As scale for 
the length y the value (5 2 ) InC0 mpr. = (fe)i as a wholly geometrical quantity, thus concentrating 
all influences of Mach number and heat transfer in the universal function 5 2 /(5 2 )«. This func- 
tion is theoretically given, when expressions for u/u§ and p/p 8 , (48) and (41), are fixed. Next 
the following term is considered: 



r d(u/u 8 ) "1 1_ / t w \ 

Ld(y/(*2)J«#« 2 \priJt 



(63) 



of eq (61) as the wall friction law for incompressible flow. For turbulent boundary layers the 
slope [d(u/us)/dy/8 2 )i] U) depends itself on the Reynolds number R S2 in such a manner that 



[ 



d(u/u 8 ) l 
(y/(82)iJw = a(H) = / r w \ 



with n=0.268 (according to Ludwieg-Tillman [18]) and 

a(y^)=0.123.10- - 678( ^ / ^ )i (65) 

as an empirical function, which expresses the influence of shape of u/u 8 . 
Therefore, we can finally write 



6 A certain confirmation for these considerations may be seen in figure 1, which shows that the velocity distribution and therefore the related 
shear stress distribution tItw too, is the same as in the incompressible case. 
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((56) 



c. The Dissipation Law 

A generalization similar to that found for the wall friction can be derived for the dissipation 
in a compressible turbulent flow with heat transfer. The energy integral condition (53) con- 
tains the term 



rdu=D (dissipation). 
Written in dimensionless form, 

PsUt PiU 2 5 J T w \UJ WphUiJi Jo T W \U8/J{8 2 ) i 



(67) 



«i.s. 



The term (T w /psUs 2 ) t is given with eq (64). The remaining integral term is found empirically 
by T. Rotta [3] and E. Truckenbrodt [4] to be 



7 ,i( u \P( H ) 



JO T w \U b / 



/>(). 1 

a{H) h ' 



(69) 



so thai il finally vi 



I) JTa(H) 0(H) ,,„ ,~| 5^ = §JJ. 



H) 5 2 
168 (5 2 )i 



(70) 



with 






0(H) 0.0052 



«* 



(71) 



as the part which represents the dissipation law for the incompressible flow. To be quite 
correct, a detailed consideration must be made on the definition of the Reynolds number 

used in eq (69). By definition the dissipation integral expresses the amount of energy dis- 
sipated by the viscosity into heat (with the intermediary of turbulence). The related Reynolds 
number should be more exactly defined with a mean viscosity jj instead of fx m because the 
energy transfer by viscosity occurs within the whole boundary layer, although mostly near 
the wall. 

Therefore eq (69) should correctly be written as 



Jo r w \uj a(H) \n/ 






(72) 



when the definition (62) for the Reynolds number is used. This correction, however, was 
proved to be in general so small that it is not essential in the scope of this approximate theory. 

d. Simplified Representation for the Wall Friction and Dissipation Law 

The expressions for the incompressible wall friction and for the dissipation law were 
given by E. Truckenbrodt [4], who evaluated some semiempirical relations found by T. Rotta [3]. 

The results of the evaluation of E. Truckenbrodt are given in his two figures 3 and 4 
which we have reproduced in the present figures 3 and 4. These have been presented in. 
older to investigate whether or not it will be possible to find a new suitable fitting for these 
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Figure 3. Wall friction law for in- 
compressible turbulent boundary 
layers. 

Thick line: evaluation of experiments of 
.J. Rotta-Truckenbrod. 

Dashed line: evaluation of experiments of 
Ludwieg-Tillmann. 

Dashed-dotted line: new approximation with 
0.2 power law. 
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Figure 4. Dissipation law for in- 
compressible turbulent boundary 
layers. 

Thick lines: evaluation of J. Rotta-Trucken- 
brodt. 

Dashed line: approximation proposed by 
Truckenbrodt. 

Dashed-dotted line: new approximation with 
0.2 power law. 



two empirical laws which permits the use (for the sake of simplification) of only one power 
of the Reynolds number instead of two (0.268 and 0.168). 

We find that if we take the power 0.2 in both cases the fitting of the experimental result 
in the most important range between R d2 =l0 3 and 10 5 will be rather satisfactory. The laws 
(66) and (70) may then be written as 



t w 0.0636 • lcr - 665 ^^' 


8. 2 


n,l/ 2_ P0.20 




(*)« 


D 0.00761 


(5 a ), 





(73) 

(74) 



This simplification prevents the universal function ^ of equation (56), from being dependent 
upon R d2 besides H, Ms and 6. Thus an important simplification in the representation of 
this function seems possible without noticable loss of accuracy. Detailed investigations on 
the reliability of this simplification are now in progress. 

This simplification will be important if the use of higher momentum equations 
(v=2, 3, 4, . . . ) is made. In this case all empirical parts of the wall friction and dissipation 
law are concentrated in the term r w /p§u 2 8 , and the integrals of type 

JoVP/ r w \uj Jo \pj bu\__T lD i>(u/u s ) J 

are dependent omV upon the shape parameter H, the Mach number Ms, and the heat transfer 
parameter d. Therefore they can be evaluated, when the family of velocity^ profiles is chosen. 
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4.3. Evaluation of the Universal Functions of the Two-Equation Method for 
Laminar and Turbulent Boundary Layers 

The evaluation of the universal functions occurring in the two-equation met hod as described 
in section 4 has already been made in the author's publications [26, 27] for the case of an insu- 
lating wall. 

The evaluation of these functions for the case with heat transfer, as suggested by the 
author, is now under way at the DVL Institute for Applied Mathematics and Mechanics, 
Freiburg/Br, by considering first the turbulent boundary layer. 

A preliminary evaluation for this case with slide rule accuracy has been available for 
about one year and was used for calculating some illustrating examples which will be discussed 
in the following section. 

This preliminary evaluation showed that the effect of heat transfer on all universal func- 
tions can be linearized in a rather good approximation, such that they may be represented 
in the form 

F,(H,M*ft) ~F v (H,M 8 0) +B'G(H } Ms). (75) 

This result will be checked after the completion of the final evaluation. 

4.4. Examples Flat Plate 

For the sake of completeness we present here again an important result obtained for 
the laminar case with insulated wall, figures 5 and 6. Pohlhausen's method based only upon 
the momentum integral condition (49) and the wall condition (13) gives H= 1.572 = const., 
for the case of dp/dx=Q at each Mach number. The consideration of the energy integral 
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Figure 5. Shape parameter H of the laminar velocity 
profile at a flat plate as a function of the Mach 
n umber Mf. 




Figure 6. Friction coefficient ratios Cf/c/ (local at the 
same Reynolds number Rs 2 ) and C/lc/ (averaged) as a 
function of the Mach number M b . 

Laminar boundary layer. 



condition (53) instead of eq (13) gives a variation of H with M* and a corresponding change 
in the behavior of the friction coefficients c s \c n (local at same values Ra 2 ) and C F /C Ft (averaged) 
with the Mach number M*. It is very interesting that the new behavior of c f /c fi (M*) is 
almost coincident with the results of an exact theory of H. Wendt [8]. This result shows 
the effectiveness of the energy integral condition. The related results for the compressible 
turbulent boundary layer with and without heat transfer are plotted in figures 7 and 8. 
The agreement with available experimental data is seen to be good. 

Figures 9 and 10 are concerned with the question of how the turbulent separation point 
in a retarded flow of 1 1 lie type 

« m Ml ' 
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Figure 7. Local friction coefficient ratio C/'c/ (at the 
same Reynolds number Rs 2 ) of a flat plate as a func- 
tion of the Mach number Ms with different values of 
the heat transfer parameter 6/r. 

Tubulent boundary layer. 
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Figure 9. Influence of the Mach number M* on 

the separation point, insulated wall (6 = 0). 

Turbulent boundary layer separation in a retarded flo\v of the type 
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Figure 8. Averaged friction coefficient ratio c/ct. 
of a flat plate as a function of the Mach number ^M b 
for insulated wall (6/r=0). 

Turbulent boundary layer. 
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is influenced by the Maeh number and by heat transfer. From figure 9 it is seen that the 

separation point (x/l) sev is almost independent of variation in the Maeh number between 
M*=0 to 2.40 (M m = to 12). Figure 10 shows that by cooling the wall in such a w.m\ 

that the wall temperature T w is kept equal to T$ ( -= 1 ), the separation is completely prevented 

even for the high value iW*=2.40 (Af 00 = 12). On the contrary, heating of the wall moves 
the separation point upstream. For the theroetical results of the figures and 10, no experi- 
mental confirmation is yet available. 

5. Conclusions 

A survey of approximation theories for compressible laminar and turbulent boundary 
layers with heat transfer at the wall and arbitrary pressure gradient in the (low direction 
shows that only a two-equation method has been used up to the present, although an infinite 
number of equations is available. A new derivation for the two infinite systems of integral 
conditions based upon the partial differential equations for the velocity and temperature 
field is given. 

An investigation of the significance of the equations with higher order shows that their 
property of describing typical boundary effects vanishes more and more with increasing number 
v. Hence, it is not surprising that the results obtained from the two-equation method with 
p=0 in both the systems (11) and (40) are good. However, the consideration of one or two 
equations more would be desirable, at least for a check of the results of the two-equation 
method in important special cases. 

Preliminary investigations have to show how the higher p-order integral terms containing 
the shear stress distribution t/t w can be evaluated in the case of compressible turbulent 
boundary layers with heat transfer. It has been found that this is possible if the power of 
the Reynolds number is chosen to be equal for both the wall friction and dissipation law. 
We have chosen 0.2. In spite of this simplification in these two empirical laws, they still fit 
the data rather well. Further work on this matter is in progress. 

This two-equation method now available with preliminary evaluated universal integral 
terms has been applied to the study of the dependence of the friction coefficient on compres- 
sibility and heat transfer. Good agreement with experimental data has been found. 
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Mechanized computation of thermodynamic tables 
at the National Bureau of Standards, J. Hilsenrath 
and J. H. Wegstein, Thermodynam. Properties of 
Fluids, London, July 10-12, 1957, Proc. Inst, of 
Mech. Engrs. (London) {1958). 

A high speed computing routine is described which computes 
the chemical composition and certain thermodynamic prop- 
erties of homogeneous gaseous systems containing up to 5 
elements. The procedure, which is applicable over a wide 
temperature range in which both dissociation and ionization 
effects are allowed for, permits the computation of mole 
fractions, compressibility factor, internal energy, enthalpy 
and entropy of the equilibrium mixtures at specified tempera- 
tures and at specified densities or pressure intervals. Auxil- 
iary routines for the IBM 704 are used to compute the 
specific heats, specific heat ratios, isentropic expansion 
coefficient and the sound velocity. A brief description is 
given of the library of thermodynamic tables that provides 
the input data for these calculations. An appendix to this 
paper gives an indexed bibliography of existing thermody- 
namic functions for the products of a wide variety of combus- 
tion reactions. 

Further identities and congruences for the coefficients 
of modular forms, M. Newman, Can. J. Math. 10, 

577 (1958). 

Let p r (n) be the coefficient of x n in 17(1— x") r . Let r be odd, 
0<r<24. Suppose p is a prime such that p>3 when(r,3) = l 
and p>2 when 3|r. Then it is shown that for all integral //. 

p r (np 2 -{-rv)— y n Pr(n)+p r - 2 p r L J j = Q) 



constants. 



(rv — ri) . 
V 



A fractional factorial experiment of the mixed 2 m 3 n 
series, W. S. Connor, Proc. U+tlx Ann. Clin. Rochester 
Soc.jor Quality Control (Feb. 1958). 

The paper is concerned with a fractional factorial design for 
three factors each of which lias two levels and two factors 
each of which has three levels. The complete factorial 
consists of 72 treatment combinations but the fractional 
factorial of only 36 treatment combinations. An example is 
worked out. 

A note on confidence intervals in regression problems 

J. Mandel, Ann. Math. Statist. No. 3, 29, 892 (Sept. 
1958). 

A method is presented for the construction of jointly valid 
confidence intervals for arbitrary real functions of the 
coefficients of a multiple linear regression. In the case of a 
single independent variable, x, the problem of constructing 
jointly valid confidence intervals for the x-values correspond- 
ing to a set of "future" measurements of the dependent 
variable is also discussed. 

Harmonic generation with ideal rectifiers, C. H. 

Page, Proc. I.R.E., No. 10, 46, 1738 (Oct. 1958). 

It is shown that the nth harmonic cannot be generated with 
an efficiency exceeding l/n 2 . Of the power converted to dc 
and harmonics, at least 75 percent is dc dissipation, and 
this cannot be reduced by an arrangement of selective 

circuits. 



Two forms of mathematical induction, A. Schach, 

Math. Mag., No. 2, 32, 83 (Nov.-Dec. 1958). 

It is shown that two principal forms of (finite) mathematical 
induction: 

(I) If P(l) and if, for all pos. int. n, P(n) implies P(/i+ 1), 

then P(ri) for all pos. int. n; and 
(II) If P(l) and if, for all pos. int. n, P(l)&P(2)&...&P(n) 
implies P{n-\-l), then Pin) for all pos. int. n 
are equivalent; and that any demonstration based on one of 
them can be recast so that it is based on the other instead. 

The evaluation of matrix inversion programs, M. 

Newman and J. Todd, J. Soc. Indust. Appl. Math., 
No. 4, 6, 466 (Dec. 1958). 

A set of thirteen test matrices typical of those occurring in 
practice were inverted on SEAC, the 704, and various other 
machines. The results were compared with the theoretical 
inverses and various error estimates were computed. 

Un nouveau critere d'univalence des transforma- 
tions dans un R n , M. A. Ostrowski, Compt. rend, 
acad. sci. 248, 348 (Jan. 1959). 

A criterion of univalence of an n-dimensional transformation, 
given in a former communication, is generalized; the possi- 
bility of further generalizations and refinements is discussed. 

On the convergence of the Rayleigh quotient itera- 
tion for the computation of the characteristic roots 
and vectors. II, A. M. Ostrowski, Arch. Rational 
Mech. Anal. No. 5, 2, 423 (1959). 

Further asymptotic properties of the Rayleigh quotient 
algorthm are derived, and numerical results for some special 
matrices are given. 

A note of the computation of x 2 , M. G. Natrella, 
Am. Statistician, No. 1, 13, 20 (Feb. 1959). 

There is a useful, but not widely-known short cut available 
for computing x 2 f° r the comparison of two observed frequency 
distributions with equal total counts. The advantages are 
that it does not require computing the expected frequencies, 
avoids problems of rounding the deviations, and requires only 
one division per class instead of one division per cell. It can 
be extended to three observed distributions, but loses its com- 
putational advantages for more than three. 

Mode expansion in the low-frequency range for 
propagation through a curved stratified atmosphere, 

H. Bremmer, J. Research NBS 63D, No. 1, 75 (1959). 

This expansion is particularly useful when considering 
ionospheric propagation at low frequencies. The complex 
problem dealing with two media, viz. a homogeneous earth 
and a surrounding stratified atmosphere, leads to intractable 
expressions. However, as the influence of the earth may 
be accounted for by an approximate boundary condition at 
the earth's surface, the problem is then reduced to that 
of the outer medium only. The coefficients of the mode 
expansion for this simplified problem will be derived while 
taking into account the earth's curvature; however, the 
latter proves to be negligible under very general conditions. 
The expansion to be derived is wanted in particular when 
studying the influence of a gradual transition in the electron 
density with height at the lower edge of the ionosphere. 
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On the theory of fading properties of a fluctuating 
signal imposed on a constant signal, H. Bremmer, 
NBS Circ. 599 (May 1959). 

This Circular is a theoretical mathematical consideration of 
the varying strength of a radio signal consisting of two 
components: a constant signal, and a fluctuating signal of 
the same frequency which has reached the receiver by a 
different path. The properties investigated for the composed 
signal are the distribution functions of both the amplitude 
and phase, as well as the average number of crossings for 
either quantity through any given level. 

The first part of this Circular deals with the idealized case, 
where the fluctuating signal is considered to be " random." 
The second part is a more practical consideration of the 
characteristics of the composite signal when the fluctuating 
signal is a " quasi-random" result of first order scattering in 
a turbulent atmosphere. In the two extreme cases, where 
the received signal is composed entirely of either the constant, 
or line-of-sight, signal, or the fluctuating, reflected signal, 
the complicated general formulas are simplified. 

This Circular will be of interest to those who deal with the 
various types of equipment and effects which pertain to the 
fading properties of this composite signal. 

The mean deviation of the Poisson distribution, 

E. L. Crow, Biometrika, pts. 3 and 4 45, 556-559 
{Dec. 1958). 

The mean deviation of^the Poisson distribution is easily 



determined to be 2mY, where m is the mean of the distri- 
bution and Y is the maximum value of the individual Poisson 
probabilities. The ratio, Rim), of the mean deviation to the 
standard deviation is shown to be a damped oscillation as a 
function of m, converging rapidly toward the limiting normal 
value V2/7T-. Furthermore Rim) takes on the normal value 
twice between every pair of successive non-negative integers. 
Thus R(m) is a poor criterion for distinguishing the Poisson 
and normal distributions. 

A property of additively closed families of distribu- 
tions, E. L. Crow, Ann. Math. Stat. No. 3, 29, 892-897 
(Sept. 1958). 

Consider a one-parameter additively closed family of univari- 
ate cumulative distribution functions F(x; X) (H. Teicher, 
Ann. Math. Stat., Vol. 25 (1954), pp. 775-778). Let three 
cumulants with orders in arithmetic progression exist and 
be nonzero. If all three orders are even, or if the first order 
is odd, it is also required that F(x; \)=0 for £<0 and F(x; A) 
> for x > 0. Consider linear combinations, with real, 
nonzero coefficients, of a finite number of independent vari- 
ables with distributions in the family. It is proved that the 
only such linear combinations whose distributions are also 
in the family are those with coefficients unity. The additively 
closed families having this property may be called strictly 
additively closed. It can be shown that (one-parameter) addi- 
tively closed stable families of distributions (normal and 
Cauchy in particular) with characteristic functions continu- 
ous in X are not strictly additively closed, while Poisson, 
generalized Poisson, binomial, and gamma families are 
strictly additively closed. 
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